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Abstract. Detrended Fluctuation Analysis (DFA), suitable for the analysis of nonstationary time 
series, is used to investigate power law in some of the Bach's pitches series. Using DFA method, 
which also is a well-established method for the detection of long-range correlations, frequency 
series of Bach's pitches have been analyzed. In this view we find same Hurts exponents in the range 
(0.7 — 0.8) in his Inventions and sinfonia. 
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INTRODUCTION 

Mathematics and music have some vague sort of affinity, but most often than not 
the supposed relationship between the two fields turns out to consist of complicated 
structure. The relationship between mathematics and musical works has been hidden 
from the listener since the old days. Thus one is forced to make use of interpretative 
techniques in order to search for them, which is problematic from a methodological 
point of view. In addition to mathematics being seen as numerical symbolism, music 
is closely linked to absolute physical entities, such as frequency and relation between 
intervals (an interval is a space between two notes). It is an illustrated fact that not just 
musical notation, but also the relationship between music and time has something to 
do with mathematics. Among great variety of complex and disordered systems most of 
music parameters such as frequency and pitch (pitch is the sound frequency of any given 
note) [11,0, HI], Amplitude or Dynamics (dynamics are the changes in volume during a 
musical piece) P], intervals (intervals are the distances between notes in the musical 
scale). Rhythm (rhythm is the structure of the placement of notes in musical time) can 
be considered as stochastic processes. Also, some authors try to cluster the music [5]. 

In this paper we characterize the complex behavior of note frequencies of a selection 
of Bach's Inventions and Sinfonias through computation of the signal parameters and 
scaling exponents. Inventions and Sinfonias are a collection of short pieces which Bach 
wrote for musical education of his young pupils. In music, an Invention is a short compo- 
sition (usually for a keyboard instrument) with two-part counterpoint, which is a broad 
organizational feature of much music, involving the simultaneous sounding of separate 
musical lines. The Inventions and Sinfonias are two and three voices music pieces re- 
spectively. These voices have independent characteristic behavior and for simplicity we 
consider only the upper voice. Because of non- stationary nature of music frequency se- 
ries, and due to finiteness of available data samples, we should apply methods which are 
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FIGURE 1. Typical up) sheet music and down) frequency series of Invention No. 6 by Bach. 



insensitive to non-stationarities, like trends. In order to separate trends from correlations 
we need to eliminate trends in our frequency data. Several methods are used effectively 
for this purpose: Detrended Fluctuation Analysis (DFA) [6], Rescaled range analysis 
(R/S) [0] and Wavelet Techniques (WT) [8]. 

We use DFA method for the analysis and elimination of trends from data sets. DFA 
method introduced by Peng et al. [6] has became a widely used technique for the 
determination of (mono-) fractal scaling properties and the detection of long-range 
correlations in noisy, non- stationary time series [| 9l. lTo|] . It has successfully been applied 
to diverse fields such as DNA sequences lS,|Tll, cardiac dynamics lfl2ll . climate [|13l1 . 
neural receptors [14], economical time series llisll etc. 

The paper is organized as follows: In section II we describe DFA methods in details 
and analyze the frequency series of the Inventions and Sinfonias. We end the paper by 
drawing conclusions. 



DFA AND ANALYSIS OF MUSIC FREQUENCY SERIES 

To implement the DFA, let us suppose we have a time series, N{i){i = 1, ...,Nmax) and 
determine the profile: y{j) = L/=i[^(0 ~ (^)]- Next we break up A'^(/) into K non- 
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FIGURE 2. a) s dependence of variation and b) The log-log plot F (s) versus s for Invention 
no. 6 frequency series. 



overlapping time intervals, !„, of equal size T where n = 0,\,...K— \ and K corresponds 
to the integer part of Nmax/ T^- In each box, we fit the integrated time series by using a 
polynomial function, ypoiii), which is called the local trend. We detrend the integrated 
time series y{i) in each box, and calculate the detrended fluctuation function: Y{i) — 
y{i) ~ypoi{i)- For a given box size s, we calculate the root mean square fluctuation: 
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The above computation is repeated for box sizes s (different scales) to provide a rela- 
tionship between F{s) and s. A power law relation between F{s) and s indicates the 
presence of scaling 

F{s)^s". (2) 

The parameter H, called Hurst exponent, represents the correlation properties of the 
signal: if H = 0.5, there is no correlation and the signal is an uncorrected signal [@]; 
if H < 0.5, the signal is anticorrelated; if H > 0.5, there are positive correlations in the 
signal. In the two latest cases, the signal can be well approximated by the fractional 
Brownian motion law [16]. Also, the auto correlation function can be characterized 
by a power law C{s) = (NkNk+s) ~ ^ '^ with 7=2 — 2H. Its power spectra can be 
characterized by S{co) ~ co^^ with frequency CO and /3 = 2H — 1. In non- stationary 
case, correlation exponent and power spectrum scaling are 7 = —2H and /3 = 2H + 1, 
respectively [[^.[itIi. 



It can be checked out that, frequency series is non- stationary. One can verify non- 
stationarity properties experimentally by measuring stability of average and variance in 
moving windows by, for example, using scale s (Fig. [2^). In Fig. [2b we plot in double- 
logarithmic scale the corresponding fluctuation function F{s) against the box size s. 



TABLE 1. Values of Hurst (//), power spectrum scaling (j3) and 
auto-correlation scaling (7) exponents for the selections of Inven- 
tions and sinfonias frequency series obtained by DFAl. 
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Invention no.l 


0.77 ±0.03 


2.54 ±0.06 


-1.54 ±0.03 


Invention no. 6 


0.72 ±0.03 


2.44 ±0.06 


-1.44 ±0.03 


Sinfonia no.l 


0.71 ±0.03 


2.42 ±0.06 


-1.42 ±0.03 


Sinfonia no. 13 


0.73 ±0.03 


2.46 ±0.06 


-1.46 ±0.03 



Using the above procedure, we obtain the following estimate for the Hurst exponent: 
H = 0.70 ±0.03. The F{s) exhibits an approximate scaling regime from 5 = 2 up to 
nearly 5 = 6 (in logarithmic scale). Since H > 0.5 it is concluded that the frequency 
series show persistence; i.e. strong correlations between consecutive increments. The 
values which derived for quantities of DFAl method for Invention no. 6 are given in 
Table [H (second line). We have calculated the Hurst exponent for other Inventions and 
Sinfonia as well, all being in the 0.7 — 0.8 range (Table [I])- 

Usually, in DFA method, deviation from a straight line in log-log plot of Eq.© occurs 
for small scales s. These deviations are intrinsic to the usual DFA method, since the 
scaling behavior is only approached asymptotically. Deviations limit the capability of 
DFA to determine the correct correlation behavior in very short records and in the regime 
of small s. The modified DFA is defined as follows |@]: 

/r^shuf/ /Ni2\l/2 1/2 

^"^°'(^) = ^(^) ;|^shuf;,;]2;i/2,n/2 (for/»i), o) 

where ([^^^"^(5)]^)^/^ denotes the usual DFA fluctuation function, defined in Eq.©, 
averaged over several configurations of shuffled data taken from original series, and 
s' ^ N/20. The improvement is very useful especially for short records or records that 
have to be split into shorter parts to eliminate problematic nonstationarities, since the 
small s regime can be included in the fitting range for the fluctuation exponent H. 

The modified DFA method indicates the correct correlation behavior also in presence 
of broadly distributed data, where the common DFA fails to distinguish long-range 
correlations from deviations caused by broad distributions. The values of Hurst exponent 
obtained by modified DFAl methods for frequency series is 0.72 ±0.03. The relative 
deviation of Hurst exponent which is obtained by modified DFAl in comparison to 
DFAl for original data is less than 4%. 



CONCLUSION 

DFA is a scaling analysis method used to quantify long-range power-law correlations 
in signals. Many physical and biological signals are 'noisy', heterogeneous and exhibit 
different types of nonstationarities, which can affect the correlation properties of these 
signals. Applying DFAl method demonstrates that the music frequency series have long 



range correlation. We calculated Hurst exponent for other Inventions and Sinfonia and 
found it to be in the 0.7 — 0.8 range. 
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